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We experimentally disentangle the contributions of different quantum paths in high-order har¬ 
monic generation (HHG) from the spectrally and spatially resolved harmonic spectra. By adjusting 
the laser intensity and focusing position, we simultaneously observe the spectrum splitting, fre¬ 
quency shift and intensity-dependent modulation of harmonic yields both for the short and long 
paths. Based on the simulations, we discriminate the physical mechanisms of the intensity-dependent 
modulation of HHG due to the quantum path interference and macroscopic interference effects. 

Moreover, it is shown that the atomic dipole phases of different quantum paths are encoded in the 
frequency shift. In turn, it enables us to retrieve the atomic dipole phases and the temporal chirps 
of different quantum paths from the measured harmonic spectra. This result gives an informative 
mapping of spatiotemporal and spectral features of quantum paths in HHG. 

PACS numbers: 42.50.Hz, 32.80.Qk, 33.80.Wz, 32.80.Wr 


I. INTRODUCTION 

Nowadays, high-order harmonic generation (HHG) 
through the interaction of intense laser pulses with 
atomic or molecular gases has been extensively inves¬ 
tigated to produce coherent extreme ultraviolet (XUV) 
radiations m and attosecond pulses mg. These ultra- 
short pulses can serve as an important tool for detecting 
the ultrafast electron dynamics inside atoms or molecules 
mm as well as inaugurating a new domain for time- 
resolved metrology and spectroscopy on attosecond time 
scale mini Eg. Due to these applications, HHG has 
been a subject of great interest in the past two decades. 

HHG process includes both the response of the single 
atom exposed to laser field and also the macroscopic co¬ 
propagation of the laser field and high harmonics. The 
physics underlying single-atom response can be well un¬ 
derstood by the three-step model M- (i) ionization, (ii) 
acceleration and (iii) recombination. Following this idea, 
a quantum theory has also been developed within the 
strong-field approximation (SFA) [15], which also can be 
incorporated into the framework of Feynman’s path in¬ 
tegral theory [16] . In this approach, the harmonic dipole 
moment can be written as the coherent sum over all dif¬ 
ferent quantum paths contributing to HHG. These quan¬ 
tum paths are a generalization of the classical electron 
trajectories. For each harmonic in the plateau region, two 
families of electron trajectories with the shortest excur¬ 
sion times (r) in the continuum, which are commonly re¬ 
ferred to the short and long quantum paths, dominantly 
contribute to the harmonic emission. The phase associ¬ 
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ated to each quantum path of the qth harmonic emission 
is given by the action of the electron accumulated in the 
laser field, which is usually called atomic dipole phase 
and can be approximated by & J q « — U p r3 « —cn? q I in 
terms of the SFA model m- Here, U p and I are the 
ponderomotive energy and laser intensity and a? is the 
phase coefficient with j = S,L representing the short and 
long quantum paths, respectively. 

The single-atom dipole phase is naturally coupled with 
the macroscopic response, which leads to different phase¬ 
matching and therefore rich spatial and spectral features 
of HHG [lS 1- Thus, in turn, the spatial and spectral fea¬ 
tures offer useful information on the HHG process [19] 20 
and the gas medium m- Moreover, each path has a spe¬ 
cific excursion time and frequency property and it leads 
to a mapping between time and frequency for each har¬ 
monic [22] which has also been used to detect the nuclear 
dynamics [23] . On the other hand, the quantum paths 
can interfere with each other, resulting in frequency split¬ 
ting f24F28] and interference fringes [291I3T] in the har¬ 
monic spectra and also laser intensity-dependent interfer¬ 
ence of harmonic yield [32ll38j . The interference fringes 
allow to retrieve even more information about electron 
motion with high precision. However, because the single¬ 
atom and macroscopic phase-matching effects are natu¬ 
rally coupled in HHG process, the mechanisms under¬ 
lying the interference and spectral features are always 
entangled and sometimes the interpretation is different 
or confusing even for the similar observations El I Mr 
33]. Experimental investigations of fully resolved HHG 
with both spatiotemporal and spectral features will shed 
more light on these issues. It is of crucial importance, 
both from a fundamental point of view, to understand the 
laser-atom interaction, as well as for applications such as 
HHG spectroscopy [39] . However, because the phase co¬ 
efficient of the long path is almost one-order of magnitude 
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higher than that of the short path, the phase-matching is 
quite different and it becomes very challenging to simul¬ 
taneously resolve the spatiotemporal and spectral fea¬ 
tures and interference fringe both for the long and short 
paths. 
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FIG. 1: (Color online) Experimental setup used for the har¬ 
monic generation and detection. 

In this work, we aim at comprehensively disentangling 
the HHG features of different quantum paths with spec¬ 
trally and spatially resolved harmonic spectra and iden¬ 
tifying the mechanisms of the intensity-dependent inter¬ 
ference of HHG. In our experiment, we observe that each 
individual high-order harmonic is spectrally split, i.e., 
shows a double-peak structure. Those two peaks are 
spatially separated in the far field and each peak also 
presents different frequency shift as increasing the laser 
intensity. By adjusting the laser intensity and focusing 
position, each peak shows different spatial profiles, ei¬ 
ther blue- or red- frequency shift and intensity-dependent 
periodic modulation. By performing the simulations 
of both the single-atom and macroscopic responses of 
HHG, we differentiate the intensity-dependent interfer¬ 
ences caused by the quantum path interference (QPI) 
and the transient phase matching. In comparison with 
simulation and experiment, the measured spectral and 
spatial features are well explained. To the best of our 
knowledge, it is the first time to observe the spectrum 
splitting, frequency shift and intensity-dependent inter¬ 
ference simultaneously for the short and long quantum 
paths in experiment. Moreover, we show that the atomic 
dipole phases of different quantum paths are encoded in 
the frequency shift. In turn, it enables us to retrieve the 
atomic dipole phases and the temporal chirps of different 
quantum paths from the measured harmonic spectra. 


II. SPATIALLY AND SPECTRALLY RESOLVED 
HIGH ORDER HARMONICS 

Our experiment is carried out using a commercial 
Tiisapphire laser system (Legend Elite-Duo, Coherent, 
Inc.), which delivers the 30-fs, 800-nm pulses at a repeti¬ 
tion rate of 1kHz with a maximum pulse energy of 10 mJ. 



FIG. 2: (Color online) (a)-(b) Measured harmonic spectra at 
the laser energies of 0.9 mJ and 3.11 mJ, respectively, (c) 
The spatially integrated HHG signals for the spectra in (a) 
(dashed line) and (b) (solid line). For the purpose of clarity, 
the solid line is multiplied by a factor of 0.2. (d) Spatially 
integrated HHG signals as a function of the laser intensity. 
Here, the gas cell is located at Z=2 mm. 


Figure [l] shows the sketch of the experimental setup. In 
brief, the incident laser beam is focused to a 2-mm long 
gas cell filled with argon with a pressure of 35 torr by a 
600-mm focal-length lens. The gas cell is mounted on a 
three axis picomotor actuator controlled by a computer. 
An argon gas jet emitted from a 500-/im diameter noz¬ 
zle has also been tried and similar phenomena (spectrum 
splitting and frequency shift) can also be observed. Of 
course, by using a gas jet, the experimental conditions 
(e.g., the gas pressure and interaction distance) are dif¬ 
ferent from that with a gas cell. A half-wave plate and 
a polarizer are placed in the light path to continuously 
control the laser energy from 0.1 mJ to 5 mJ, and an iris 
diaphragm is used to adjust the laser beam size (around 
9-mm diameter in the experiment). A 500-nm thick alu¬ 
minum foil is placed in front of the spectrometer to block 
the driving pulse. The generated harmonic spectrum is 
detected by a home-made flat-field soft x-ray spectrom¬ 
eter consisting of a flat-field grating (1200 grooves/mm) 
and a slit with a width of about 0.1 mm and hight of 15 
mm. High order harmonics are dispersed by the grating 
and imaged onto the Micro-channel plate (MCP) fitted 
with a phosphor screen. The image on the screen is read 
out by a charge-coupled device (CCD) camera. 
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FIG. 3: (Color online) (a) Measured harmonic spectrum with the gas cell located at Z=—3 mm, (b) spatially integrated 
harmonic spectra as a function of the laser energy, and (c) calculated spatially resolved harmonic spectrum to compare with 
(a), (d)-(f) and (g)-(i) are the same to (a)-(c), but for the gas cell located at Z=0 mm and Z=3 mm, respectively. 


Figures Ha) and [2jb) display the spatially resolved 
harmonic spectra obtained at the laser energies of 0.9 
mJ and 3.11 mJ respectively, and their corresponding 
spatially integrated HHG signals are presented as the 
dashed (0.9 mJ) and solid (3.11 mJ) lines in Fig. 2(c). In 
our experiment, we have estimated the laser intensity by 
measuring the laser power and the focus size. The laser 
power can be directly read out from a power meter. The 
focus size is measured by using a CCD beam profiler that 
mounted on a translation stage in the direction of the 
beam propagation, and finally the measured focal spot is 
demonstrated to have a diameter about 200 fi m. Due to 
the dispersion of the optical elements (e.g., a 3-mm thick 
window plate in the front of the vacuum chamber and a 
3-mm thick half-wave plate), the 30-fs fundamental laser 
pulse delivered by the laser system will be broadened 
to about 36 fs. Finally, the laser intensity is estimated 
about 2.75 x 10 14 W/cm 2 and 0.8 x 10 14 W/cm 2 at the 
laser energy of 3.3 mJ and 0.9 mJ, respectively. How¬ 
ever, with this method, the uncertainty of the estimated 
laser intensity may be up to 20%. The laser intensity 
can also be determined from the cutoff of the measured 
harmonic spectra. However, due to the limitation of the 
dynamic range of CCD, the cutoff position is difficult to 
identify from the generated harmonic spectra. The gas 
cell is placed at Z= 2 mm. Here Z< 0 and > 0 means 
the upstream and downstream of the laser focus, respec¬ 
tively. One can clearly see several remarks from Fig. [2| 


First, each individual harmonic is split into two branches 
both in the spatial and the spectral domains. The spatial 
and spectral separations at high laser intensity (3.11 mJ, 
- 2.75 x 10 14 W/cm 2 ) are much larger than those at a low 
laser intensity (0.9 mJ, ~ 0.8 x 10 14 W/cm 2 ). Second, 
the spatial divergence of left branch is much broader (±6 
mrad) and depends more sensitively on the laser inten¬ 
sity than that of the right one (±2 mrad). This suggests 
us to assign the left branch as the long path contribution 
and the right one as the short path contribution, since 
the long path requires a lager divergence angle for phase 
matching than the short path [40] . Third, as shown in 
Fig. Hd) , the left branch of each individual harmonic 
(long path) is clearly blue-shifted, while the right branch 
(short path) is red-shifted in the frequency domain as 
the laser intensity increases. For simple description, in 
the following, we abbreviate the blue and red shift of the 
short (long) path as SB (LB) and SR (LR) shift, respec¬ 
tively. 

In Fig. [3] we have presented the harmonic emission 
at different focusing geometries. Figure |3ja) shows the 
spatially resolved harmonic spectrum generated with the 
laser intensity of 2.3 x 10 14 W/cm 2 (2.6 mJ) at Z=—3 
mm. Figure [3^b) shows the spatially integrated harmonic 
spectra as a function of laser intensity. All these indi¬ 
vidual harmonics are spatially and spectrally separated. 
However, the spatial and spectral profiles in Figs, 
and (b) are quite different from Fig. [ 2 ] One can see 
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FIG. 4: (Color online) (a) Temporal profile of the laser in¬ 
tensity. (b) Instantaneous frequencies of H15 for the short 
(dashed line) and long (solid line) paths. (c) Instanta¬ 
neous frequencies of HI5 as a function of the laser intensity, 
(d) Comparison between experiment and simulation for the 
intensity-dependent 15th harmonic emission. 


that, at Z=—3 mm [Figs. !a), (b)], the left branch is 
distributed in a narrower space range (±2 mrad) and 
blue-shifted in the frequency domain. While the right 
branch presents a broader spatial distribution within 
—5^3 mrad and a red-shift. In other words, the left and 
right branches in Figs. [3|a) and (b) show the SB and LR 
shift. In contrast, those in Fig. [2] show the LB and SR 
shift. Distinct and interesting spatial and spectral fea¬ 
tures also can be observed at Z=0 mm [Figs. [3] (d) and 
(e)] and Z=3 mm [Figs. [3j[g) and (h)]. To be specific, 
for Z=0 mm, both the left and right branches are blue- 
shifted. However, the left branch shows a large spatial 
divergence (i.e., LB) and the right branch shows a small 
spatial divergence (i.e., SB). In contrast, for Z=3 mm, 
each of the two branches has a narrow spatial divergence 
(±2 mrad). However, the left branch is blue-shifted (i.e., 
SB) and the right branch is red-shifted (i.e., SR). In ad¬ 
dition, we have also calculated far-field spatial profiles 
of harmonics at Z=—3 mm, Z=0 mm, and Z=3 mm to 
compare with the experimental results in Figs. [3] (a), (dV, 
and (g). Corresponding results are presented in Figs, p] 
(c), (f), and (i), respectively. One can see that, the calcu¬ 
lated results agree reasonablely with the measurements: 
each harmonic is split into two branches, and the spatial 
distributions of the separated two branches (attributed 
to short or long path) are in reasonable agreement with 
the experiments. Some differences in details are likely 
due to the uncertainties of the experimental conditions, 
e.g., the laser intensity and focus size. In short, by ad¬ 
justing the laser intensity and focusing position, we have 
simultaneously observed rich spatio-spectral features of 
the spectrum splitting and frequency shift both for the 
short and long paths, which have not been observed so 
far. 


The observed spectrum splitting and frequency shift 
can be explained by considering the macroscopic phase¬ 
matching of HHG j4QH42| . The phase mismatch between 
the driving laser field and the emitted high order harmon¬ 
ics can be expressed by A k = A k q + A k g + A k a + A k e . 
Here A k q is the phase mismatch due to the single-atom 
dipole phase and can be written as —a q VI. A k g is the 
geometrical wave vector mismatch due to focusing. A k a 
and A k e are the wave vector mismatches between high 
harmonics and fundamental field due to the dispersion in 
the neutral medium and free electrons. Due to the vari¬ 
ation of laser intensity with time, A k q and therefore A k 
will change during the laser pulse. This leads to transient 
phase matching of different quantum paths. The varia¬ 
tion of the instantaneous frequency of each harmonic can 
be expressed by omum 

u q {t)=qu o + (1) 

In Fig. gb), we show the instantaneous frequencies of 
the short (dashed line) and long (solid line) paths for the 
15th harmonic (HI5) as an example. One can see that at 
the uphill of the laser pulse, both the short and long paths 
present a blue shift, i.e., the SB and LB shift. While at 
the downhill, red shift is presented for these two paths, 
i.e., the SR and LR shift. Since the phase coefficient 
a q a q , the frequency shift of the short path is much 
weaker than the long one, thus is more difficult to be ob¬ 
served in experiments. During the laser pulse, the phase 
matching can be transiently optimized at certain inten¬ 
sities I m due to the intensity dependence of the atomic 
dipole phase. By inserting a Gaussian intensity envelope 
(i.e., I(t) = Ioexp(— 41n2t 2 /r 2 ), where Jo and r denote 
the peak intensity and the pulse duration full width at 
half maximum) into Equation (1), one can then calculate 
the Jo-dependent instantaneous harmonic frequency at a 
certain phase matching intensity J m by [33l 05]: 

u q (J 0 ) = quo ± oP q fAFE yintlo/I m ). (2) 

r 

Figure |4][c) shows the calculated intensity dependence 
of the instantaneous frequency of H15 with J m =2.0 x 
10 14 W/cm 2 . With the increase of the laser inten¬ 
sity, HI5 is gradually separated into four branches, as 
marked as the 1(LB), 2(SB), 3(SR), and 4(LR), respec¬ 
tively. When compared with the experimental result at 
Z=2 mm, we find that evolutions of the measured two 
branches match well with the branch 1 (LB, left) calcu¬ 
lated with J m =0.23 x 10 14 W/cm 2 (0.26 mJ) and branch 
3 (SR, right) with J m =1.34 x 10 14 W/cm 2 (1.52 mJ), 
respectively. In our experiment, we have observed vari¬ 
ous combinations of the four branches at different focus¬ 
ing positions. However, it is very challenging to capture 
these four branches simultaneously, because the phase co¬ 
efficient of the long path is almost one-order of magnitude 
higher than that of the short path. 

Tracing the four branches requires strict experimental 
conditions to satisfy the phase matching of both long and 
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FIG. 5: (Color online) The harmonic spectrum measured at 
the intensity of 2.9 x 10 14 W/cm 2 (3.3 mJ) with the gas cell 
placed at Z=—6 mm. 

short paths and also large frequency splitting. To this 
end, we have tried to observe the harmonic spectrum at 
Z=—6 mm (to ensure the coexist of the short and long 
paths) with a high intensity about 2.9 x 10 14 W/cm 2 (3.3 
mJ, to make the two paths both well separated). Under 
this condition, the ionization probability is so high that 
the harmonic yield decreases rapidly due to the disper¬ 
sion induced by free electrons and the absorption in the 
gas medium. Figure [5] shows the spatially resolved har¬ 
monic spectrum. One can see that the 15th and 17th har¬ 
monics have been split into 3 (or 4) branches in the spa¬ 
tial and spectral domains. The outer most two branches 
exhibit relatively broader spatial profiles compared to the 
inner one, which are mainly due to the separation of the 
long path, corresponding to the branches 1 and 4, re¬ 
spectively. While the inner one is the coalition of the 
branches 2 and 3 from the short path. 



Wavelength (nm) 


FIG. 6: (Color online) Spatially integrated HHG signals as a 
function of the gas pressure. Here, the gas cell is located at 
Z=2 mm, and the laser intensity is fixed at 1.77 x 10 14 W/cm 2 
(2 mJ). 

In our experiment, the phase matching conditions of 
different quantum paths are changed by adjusting the 
laser intensity and focusing position, which mainly affects 
the values of A k q and A k g and therefore changes the A k. 
When the laser intensity is high enough, large numbers 
of electrons are ionized. The mismatching due to the 
dispersion of free electrons (A k e ) as well as the depletion 
of the ground state becomes serious and cannot be ne¬ 


glected, which will consequently lead to a rapid decrease 
in the harmonic yield (as mentioned in Fig. [5] above). 
However, in this case, the spectrum splitting can still be 
observed from the measured harmonic signals (as shown 
in Fig. [5]). For the dispersion in the neutral medium 
A k a , it depends on the gas pressure. In Fig. [6j we show 
the measured harmonic spectra as a function of the gas 
pressure. Here, the gas cell is located at Z=2 mm, and 
the laser intensity is fixed at 1.77 x 10 14 W/cm 2 (2 mJ). 
With the gas pressure changing from 10 to 50 torr, the 
spectrum splitting of each harmonic can always be ob¬ 
served, but the strengths of the separated two branches 
increase gradually. These results indicate that the ion¬ 
ization and dispersion of the gas medium (A k e and A k a ) 
could have an effect on the measured harmonic strength, 
but will not affect the spectrum splitting of the generated 
harmonics. 

On the other hand, it should be noted that the sepa¬ 
rated branches in harmonic spectra show obvious asym¬ 
metry. For instance, in Fig. [4] (d), the harmonic emission 
from 1LB is stronger than that from 3SR, especially in 
the intensity range from 0.88 x 10 14 W/cm 2 (1 mJ) to 
1.77 x 10 14 W/cm 2 (2 mJ). And in Fig. [HJ the branch 4 
(LR) is much brighter than branch 1 (LB) and dominate 
the harmonic emission. To clarify this asymmetry of the 
harmonic splitting, we next investigate the phase match¬ 
ing of these separated branches. The phase mismatch A k 
for both the short and long paths can be calculated with 
the theory in [40] . In the following, we take the 15-th har¬ 
monic as an example. Under our experimental condition, 
the refractive indexes of the laser field and harmonics in 
the gas are very close and therefore the phase mismatch¬ 
ing due to the dispersion in the neutral medium (A k a ) is 
very small and can be neglected. In the case of Fig. [4](d), 
the dominant HHG contributions are from the long path 
with a blue shift (1LB) and the short path with a red shift 
(3SR). Then we analysis the phase matching of the long 
(short) path at the uphill (downhill) of the laser pulse. 
The laser intensity here is chosen as 1.5 x 10 14 W/cm 2 to 
ensure the asymmetry is distinct. Under this intensity, 
the ionization probability of electrons is below 2%, the 
electronic dispersion (A k e ) is small and can also be ne¬ 
glected. We have calculated the phase mismatching A k q 
and Ak g for the long path at t = —10 fs (uphill) and 
the short path at t = 10 fs (downhill). The result shows 
that, for the long path, Ak q =— 264.6/m, A/^=756.1/m, 
and the total phase mismatching A&=491.5/m. While 
for the short path, Ak q =— 16.5/m, A/^=756.1/m, and 
A/c=739.6/m. One can clearly see that for the short 
and long paths, A k g is the same, but A k q is quite dif¬ 
ferent. Finallythe total phase mismatching of the long 
path is smaller than that of the short one. As a conse¬ 
quence, the harmonic emission from the branch 1LB is 
stronger than 3SR. This result indicates that the asym¬ 
metry of the two branches in Fig. [4] (d) is mainly due 
to the different single-atom dipole phases of the short 
and long paths. Next we consider the phase matching 
in a more intense laser field, e.g., 2.9 x 10 14 W/cm 2 as 
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in Fig. [5] In this case, the ionization probability is 
up to 50%. The electronic dispersion term (A k e ) can 
no longer be ignored. We have calculated the phase 
mismatching A k q , A k g , and A k e for the long path at 
t=—10 fs (uphill, for branch 1LB) and t=10 fs (down¬ 
hill, for branch 4LR). For branch 1LB, A/c g =1310/m, 
Afc p =698.7/m, and Ak e =— 6.5/m. While for 4LR, 
A/c g =1310/m, AA^=698.7/m, and Ak e =— 1050/m. One 
can see that both A k q and A k g are equal for these two 
branches. But the electronic dispersion (A k e ) becomes 
more serious at t=10 fs. This is due to the larger density 
of free electrons at the downhill of the laser pulse. The 
total phase mismatching A k of branch 1LB is 2002.2/m, 
whereas it is only 958.7/m for branch 4LR. In other 
words, the electronic dispersion A k e can compensate the 
mismatching due to the dipole phase and geometrical fo¬ 
cusing, the total mismatching of the long path is smaller 
at the downhill of the laser pulse. Therefore, the branch 
4LR is brighter than the branch 1LB and dominate the 
harmonic emission. 


III. INTENSITY-DEPENDENT 
INTERFERENCE OF HHG YIELD 
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FIG. 7: (Color online) Spatially integrated harmonic signals 
of H15 and H21 generated at Z=2 mm with (a) no spatial 
filtering, (b) off-axis spatial filtering (from 2 to 6 mrad), and 
(c) on-axis spatial filtering (from —2 to 2 mrad). For a better 
comparison, the dashed line is multiplied by a factor of 3. 

Apart from the spectrum splitting and the frequency 
shift discussed above, one can also observe pronounced 


periodic modulations of HHG yields from Fig. HP). In 
detail, we have presented in Fig. [ 7 ] the intensity de¬ 
pendence of HHG yields for HI5 (solid lines) and H21 
(dashed lines) with the gas cell placed at Z=2 mm. Fig¬ 
ure [TJa) shows the result obtained without any spatial 
filtering, where both the short and long paths are in¬ 
cluded. As shown in this figure, these two harmonics 
show the same modulation periodicity of about 0.3 ^ 
0.4 x 10 14 W/cm 2 . In Fig. [7|b), the off-axis spatial fil¬ 
tering (from 2 to 6 mrad) is performed such that the long 
path is selected. The obtained modulations are basically 
the same as in Fig. [7|a), but the modulation amplitude 
is reduced. Similar results are also found in Fig. [TJc), 
where a spatial filter on axis (from —2 to 2 mrad) is used 
such that the short path dominates the harmonic emis¬ 
sion. These results indicate that the modulation of HHG 
yields in our experiment depends insensitively on the har¬ 
monic order as well as the dominant quantum path. 



FIG. 8: (Color online) Calculated HHG yields with the single¬ 
atom response. Solid and dashed lines are for HI5 and H21 
with both the short and long paths contributions. The dotted 
and dash-dotted lines are for the 15th harmonic with only the 
long or short path contribution. Here, S and L represent the 
short and long paths. 

Similar intensity-dependent modulations of HHG 
yields have been previously reported [32l 155 ] . However, 
the interpretation is different and still confusing. At the 
single-atom level, the modulation has been attributed to 
the QPI between the short and long paths due to their 
different intensity-dependent dipole phases [32] . While 
on the macroscopic aspect, the transient phase match¬ 
ing due to the variation of Ak with time can also lead 
to oscillation of HHG yield, which is also called tem¬ 
poral Marker fringes mmm- To disentangle these 
two physical mechanisms, we perform the simulations of 
harmonic emission based on both the single-atom and 
macroscopic phase-matching responses. Figure [8] shows 
the calculated HHG yields with the single-atom response 
[32] . It is obvious that when the short and long paths co¬ 
exist, QPI occurs and induces distinct modulations in the 
intensity-dependent HHG yields for both HI5 (solid line) 
and H21 (dashed line). Whereas, the modulation peri- 
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FIG. 9: (Color online) (a)-(c) Time-frequency spectrogram of the harmonic spectra generated at 2.5 x 10 14 W/cm 2 by putting 
the gas cell at the focus (Z=0 mm), the downstream (Z=2.5 mm), and the upstream (Z=—2.5 mm), (d)-(f) Intensity dependence 
of HHG yields calculated for different harmonic orders. 


ods for these two harmonics are slightly different (about 
0.25 x 10 14 W/cm 2 for H15 and 0.3 x 10 14 W/cm 2 for 
H21). This is caused by the different Aa q (difference 
of the two phase coefficients, i.e., — a\ q ) of H15 and 

H21 and the modulation period is approximately equal 
to 2n/Aa q . In Fig. [81 the intensity-dependent HHG 
yield of H15 with only the long (short) path contribution 
is also presented as the dotted (dash-dotted) line. With 
only one quantum path contribution, QPI no longer oc¬ 
curs and therefore the periodic modulation disappears. 
This simulation shows that the intensity-dependent in¬ 
terference caused by QPI demands the coexistence of the 
short and long paths and the interference fringe depends 
on the harmonic order. 

The macroscopic response in the gas medium is de¬ 
scribed by the copropagation of the laser and high har¬ 
monic fields, which is simulated by numerically solving 
the Maxwell wave equations as in [47, 48] . Similar to our 
experiment, the gas medium position is changed in the 
simulations to satisfy the phase-matching of the short, 
long and both quantum paths. Figs. c) display 

the time-frequency spectrogram calculated by the Ga¬ 
bor transform of the generated harmonic spectra 


at 2.5 x 10 14 W/cm 2 by putting the gas cell at the fo¬ 
cus (Z=0 mm), the downstream (Z=2.5 mm), and the 
upstream (Z=—2.5 mm), respectively. One can clearly 
see that, by putting the gas cell at the focus, both the 
short and long paths [corresponding to the lower and up¬ 
per branches of each peak in Fig. |9][a)] exist in each 
half optical cycle. While at the downstream, only the 
short (lower branch) path is dominant. While at the up¬ 
stream, the long path (upper branch) is well selected. 
We have calculated the HHG yields (spatially integrated 
in the far field) for different harmonic orders at the 
above three focusing positions. Corresponding results 
are shown in Figs. mm respectively. As shown 
in Fig. !<*), due to the QPI of the short and long 
paths, one can clearly see the intensity-dependent mod¬ 
ulations of HHG yields. The modulation period depends 
on the harmonic order. It is about 0.25 x 10 14 W/cm 2 for 
H15, 0.3 x 10 14 W/cm 2 for H21, and 0.4 x 10 14 W/cm 2 
for H27. In Fig. M[e), the harmonic emission is dom¬ 
inated by the short path. The calculated HHG yields 
first increase rapidly as the laser intensity changes from 
0.5 x 10 14 W/cm 2 to 2 x 10 14 W/cm 2 . In this range, the 
intensity modulation is unapparent. However, when the 





























laser intensity is over 2 x 10 14 W/cm 2 , obvious modu¬ 
lation structures emerge (see the illustration). For dif¬ 
ferent harmonic orders (here, H15, HI7, and H19 are se¬ 
lected as examples), almost the same modulation period 
of about 0.3 ~ 0.4 x 10 14 W/cm 2 is shown. Note that 
the QPI does not work in this case and the modulation 
is mainly caused by the transient phase matching of the 
short path. Similar results can also be found in Fig. [9^f), 
where the long path contribution is dominant. Here, it 
is worth mentioning that in Fig. [9^e), one can also see 
the modulations in the intensity from 0.5 x 10 14 W/cm 2 
to 2 x 10 14 W/cm 2 . This is mainly due to the imper¬ 
fect selection of the long path. These results show that 
the intensity-dependent interference fringe caused by the 
macroscopic phase-matching depends less sensitively on 
the harmonic order. This provides a criteria to differ¬ 
entiate the interference caused by QPI and macroscopic 
phase-matching effects. Recall that even though short 
and long paths coexist and the modulation depth is in¬ 
deed smaller as comparing Fig. J7^b) with[7ja), the modu¬ 
lation period depends insensitively on the harmonic order 
and dominant quantum path. Therefore, the observed 
intensity-dependent modulations in our experiment can 
be mainly attributed to the macroscopic phase-matching 
of each individual quantum path, which probably is dom¬ 
inant over the QPI. 



FIG. 10: (Color online) (a) Experimentally retrieved dipole 
phase coefficients of the short and long paths, (b) Calculated 
traveling times of the short and long paths. 


IV. RETRIEVING THE ATOMIC DIPOLE 
PHASE AND TEMPORAL CHIRP 

In section II, we have shown that the spectrum split¬ 
ting and frequency shift are mainly induced by the tran¬ 
sient phase-matching due to the variation of the atomic 
dipole phases of different quantum paths in time. In 
turn, we can retrieve the atomic dipole phase via the 
intensity-dependent frequency shift. Moreover, the spa¬ 
tially and spectrally resolved HHG provides the informa¬ 
tion for both the short and long paths simultaneously. 
According to Eq. (2), we know that when the laser inten¬ 
sity Jo reaches the phase matching intensity J m , the fre¬ 
quency shift of the qth harmonic approaches to zero. This 
allows us to directly read out the J m (where the frequency 


shift is minimum) for each harmonic from the measured 
intensity-dependent harmonic spectra. As shown in Fig. 
Id). for the short paths of H13—H19 (right branches), 
the transient phase-matching intensities can be directly 
read out at the minimum frequency shift, which are 
1.4 x 10 14 W/cm 2 (1.59 mJ), 1.34 x 10 14 W/cm 2 (1.52 
mJ), 1.3 x 10 14 W/cm 2 (1.47 mJ), and 1.25 x 10 14 W/cm 2 
(1.41 mJ), respectively. Then one can retrieve the phase 
coefficient for each harmonic (<af 3 = 4.6 x 10“ 14 cm 2 /W, 
<af 5 = 7.2 x 10 -14 cm 2 /W, af 7 = 10 x 10 -14 cm 2 /W, 
and af 9 = 12 x 10 -14 cm 2 /W) by fitting the intensity- 
dependent red-shift value according to Eq. (2). In prin¬ 
ciple, the phase coefficient a q of the long quantum path 
can be retrieved in the same way. However, the transient 
phase-matching intensity of the long path is about one- 
order of magnitude smaller than that of the short one, 
i.e., < 0.4 x 10 14 W/cm 2 , where the harmonic signal and 
cutoff energy are very low. Hence it makes it difficult 
to directly obtain the J m as the short path. Then, we 
retrieve the phase coefficient by seeking the best fit value 
by setting a q in the range of [15, 25] xlO -14 cm 2 /W 
and I m in the range of [0.1, 0.4] xlO 14 W/cm 2 . Fi¬ 
nally, the optimal values of I m of the long paths for 
H13—H19 are determined as 0.18 x 10 14 W/cm 2 (0.2 mJ), 
0.23 x 10 14 W/cm 2 (0.26 mJ), 0.27 x 10 14 W/cm 2 (0.3 
mJ), and 0.31 x 10 14 W/cm 2 (0.35 mJ). Corresponding 
retrieved phase coefficients are <af 3 = 24 x 10“ 14 cm 2 /W, 
af 5 = 22.6 x 10- 14 cm 2 /W, ot\ 7 = 20 x HT 14 cm 2 /W, 
and <afg = 18.5 x 10 -14 cm 2 /W, respectively. Figure 
[To|[a) plots the retrieved a q of the short (circles) and 
long (squares) paths as a function of harmonic orders. 
Clearly, a q increases as harmonic order increases for the 
short path. In contrast, a q decreases for the long path, 
which is in agreement with those extrapolated from the 
chirp measurements [33EE] and theoretical calculations 
[52l [53]. According to the SFA model, the phase coeffi¬ 
cients are directly linked to the excursion times of quan¬ 
tum paths in terms of & « — U p r^ ~ —a? q I. For com¬ 
parison, the calculated excursion times of the short and 
long paths are also presented in Fig. The experi¬ 

mental curves show very similar trends to the theoretical 
calculations as the harmonic order increases: the short 
path shows a positive temporal chirp while the long path 
shows a negative one. 


V. CONCLUSION 

In conclusion, we have simultaneously observed 
the spectrum splitting, frequency shift and intensity- 
dependent interference both for the short and long paths 
in experiment. It is shown that each individual high- 
order harmonic is gradually split into two branches in 
the spatial and spectral domains as the laser intensity 
increases. By adjusting the focusing position and laser 
intensity, each branch shows distinct spatial profiles, ei¬ 
ther blue- or red-frequency shift. The HHG yield also 
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shows an intensity-dependent interference fringe. The 
simulations indicate that the interference fringe due to 
QPI depends on the harmonic order. In contrast, the 
interference fringe due to the transient phase matching 
of individual quantum path depends less sensitively on 
the harmonic order. This provides a robust criteria to 
differentiate these two interference mechanisms. More¬ 
over, the atomic dipole phase is encoded in the intensity- 
dependent frequency shift, from which we retrieve the 
atomic dipole phases and reveal the temporal chirps both 
for the short and long paths. 

The heart underlying the rich spatiotemporal and spec¬ 
tral features is the atomic dipole phase = —a q I(t) 
associated with different quantum paths. The differ¬ 
ent phase coefficients lead to spatially resolved harmonic 
signals for short and long paths. The evolution of the 
intensity-dependent dipole phase in the time domain 
leads to transient phase matching, which results in spec¬ 
tral splitting, intensity-dependent frequency shift and in¬ 
terference of HHG yields. The frequency shift enables us 
to measure the phase coefficient a q , which depends on the 
harmonic order. In turn the order-dependent phase coef¬ 
ficients reveal the positive and negative temporal chirps 
of the short and long paths, i.e., the time-frequency prop¬ 
erty of HHG in the sub-cycle of the laser field. 

Previous attosecond spectroscopy techniques are 


mainly concentrated on the spectral features of harmonic 
spectra or the temporal features of the generated at¬ 
tosecond pulses. On the other side, as coherent lights, 
high harmonics contain other intrinsic features, such as 
the spatial profiles, interference fringes, etc.. A few re¬ 
cent investigations have proposed and demonstrated that 
these new features can shed new light on the molecular 
internal dynamics [36] and in situ measurement of iso¬ 
lated attosecond pulses in space and time |39]. Our work 
provides an informative mapping of the spatiotemporal 
and spectral features of quantum paths in HHG. These 
observed HHG features can serve as promising tools for 
attosecond spectroscopy in both spectral and spatiotem¬ 
poral dimensions. 
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